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ABSTRACT 

We study the relation between the chemical composition and the type of dust present in a 
group of 20 Galactic planetary nebulae (PNe) that have high quality optical and infrared spec¬ 
tra. The optical spectra are used, together with the best available ionization correction factors, 
to calculate the abundances of Ar, C, Cl, He, N, Ne, and O relative to H. The infrared spectra 
are used to classify the PNe in two groups depending on whether the observed dust features 
are representative of oxygen-rich or carbon-rich environments. The sample contains one ob¬ 
ject from the halo, eight from the bulge, and eleven from the local disc. We compare their 
chemical abundances with nucleosynthesis model predictions and with the ones obtained in 
seven Galactic H II regions of the solar neighbourhood. 

We find evidence of O enrichment (by ~ 0.3 dex) in all but one of the PNe with carbon- 
rich dust (CRD). Our analysis shows that Ar, and especially Cl, are the best metallicity in¬ 
dicators of the progenitors of PNe. There is a tight correlation between the abundances of 
Ar and Cl in all the objects, in agreement with a lockstep evolution of both elements. The 
range of metallicities implied by the Cl abundances covers one order of magnitude and we 
find significant differences in the initial masses and metallicities of the PNe with CRD and 
oxygen-rich dust (ORD). The PNe with CRD tend to have intermediate masses and low metal¬ 
licities, whereas most of the PNe with ORD show higher enrichments in N and He, suggesting 
that they had high-mass progenitors. 

Key words: ISM: abundances - planetary nebulae: general - H II regions - stars: abundances 
- nuclear reactions, nucleosynthesis, abundances. 


1 INTRODUCTION 

Planetary nebulae (PNe) are the final products of many stars with 
masses below ~ 8 Mq. This evolutionary stage occurs after the 
asymptotic giant branch (AGB) phase, if the star reaches a tempera¬ 
ture high enough to ionize the circumstellar gas (the AGB external 
envelope ejected through a stellar wind) before it dilutes into the 
interstellar medium. The ionized gas contains valuable information 
on the nucleosynthesis processes occurring inside the AGB stars 
(because some elements are produced in the interior of these stars 
and then carried to the surface) and on the chemical composition 
of the environment where the star was bom (since other elements 
remain unchanged during the life of the star). The chemical abun¬ 
dances in PNe can be compared with the ones computed using the 
same techniques in H II regions in order to obtain information not 
only about the efficiency of stellar nucleosynthesis and dredge-up 
processes but also about the development of chemical evolution in 
different galaxies. 

The oxygen abundance has been traditionally used as a proxy 
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for the metallicity in ionized nebulae because oxygen is the ele¬ 
ment for which more reliable abundances can be obtained. Bright 
emission lines from two of its ionization states, and can be 
easily measured in optical spectra, and the correction for the con¬ 
tribution of higher ionization states is large only for PNe of very 
high excitation. However, AGB stars can modify the oxygen abun¬ 
dance via two mechanisms. The third dredge-up (TDU) is a mixing 
event that transports to the stellar surface material created by He 
burning and the .s-process. The materi al is heavily enriched in car - 
bon but also has some oxygen (see e.g. lKarakas & Lattanziol2()I^ . 
On the other hand, hot bottom burning (HBB) occurs in the most 
massive progenitors and consists in the penetration of the bottom of 
the convective envelope into a region where the temperature is suf- 
fici ent to activate the CNO-cy cle, which can destroy oxygen (see 
e.g. lKarakas & LattanziollOldh . The products are then conveyed to 
the surface. 

The efficiencies of the TDU and HBB depend on the initial 
mass of the star and on its metallicity, and the predictions dif¬ 
fer from model to model. For example, standard nucleo synthesis 
model s with no extra mixing processes, such as those of lKarak^ 
koiol) . do not predict a significant production of oxygen at solar 
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or half solar metallicity. However, models that i nclude d i ffusive 
con vective overshooting, like those computed by iMarigol l l200lh 
and IPignatari et al]ll2013ll . predict a significant production of oxy¬ 
gen, even at s olar metallicit i es. Th is mixing mechanism was first 
introduced by iHerwig et alj l ll997h and causes the convective en¬ 
velope to penetrate more deeply into the star and, in combination 
with an efficient third dredge up, leads to an increase of the oxy¬ 
gen abundance for stars between 1.5 Mq and 3 Mq. There are other 
mechanisms that can produce extra mixing, s uch as rotation, mag¬ 
netic fields, and thermohaline mixing (see e.g. lKarakas & Lattanziol 
l2014h . These mechanisms are not well understood, but reliable de¬ 
terminations of the chemical composition of ionized nebulae might 
help us to constrain their efficiencies. Observationally, claims of 
oxygen production in l ow mass stars have been restricted to low 
metal li city PNe (see e.g.lFe guignot et al.l20()3: iLeisv & Dennefel^ 
l2006h . iRodrfguez & Delgado-Ingladal j201lh found an overabun¬ 
dance of oxygen in nearby PNe of near-solar metallicity when they 
compared them with local H II regions, but a detailed analysis of 
the clues provided by other elements is required in order to ascer¬ 
tain the c ause of this ove rabu ndance. On the other h and, the AGB 
models of lKarak^ l l2010tl and IPignatari et '^l l2013l) do not predict 
a significant destr uction of oxygen for metallicities above Z ^ 0.01, 
but the models of iKarakall ( 1201 Oh with Z 0.008 show that 0/H 
can decrease by ~ 0.1-0.2 dex for stars with masses M > 5 Mg, 
and this effect could also be m ore important at lower metallicities 
(see e.g. lStasinska et alj201(]h . 

Other elements whose abundances we can derive and use as 
proxies for nebular metallicities are sulfur, neon, chlorine and ar¬ 
gon. However, these elements require corrections for the contri¬ 
bution of unobserved ions to their total abundances, the so-called 
ionization correction factors (ICFs). Up until recently, the avail¬ 
able ICFs were based on the similarities of ionization potentials 
of the ions involved or on the r esults provided by a small n umber 
of photoionization models (e.g. lKingsburgh & Barlowlll994h . This 
proble m is now alleviated by the new and more reliable ICFs calcu¬ 
lated bv IPelgado-Inglada, Morisset. & Stasinskal ( l2014h . which are 
based on large grids of models. These new ICFs can be used to de¬ 
rive more reliable estimates of the chemical abundances, including 
those elements that are expected to be modified during the evolu¬ 
tion of the stellar progenitors of PNe, namely, helium, nitrogen, and 
carbon. 

A very useful abundance ratio in the study of PNe is C/0. Its 
value is related to the initial mass and metallicity of the progenitor 
star since these parameters determine the efficiencies of the nucle¬ 
osynthesis processes that change the relative abundances of C and 
O. Stars are bom with C/0 < 1 (C/O - 0 .5 for solar metallicity; 
lAllende Prieto. Lambert. & AsDlundl2002h . The lowest mass AGB 
stars (with initial masses below 1.5 Mq), do not go through the 
TDU or HBB, and thus keep their C/O values below 1. They are 
thus defined as 0-rich stars. Intermediate mass stars, with masses 
between ~ 1.5-2.0 Mg and ~ 4-5 Mg (the mass range depends 
on metallicity and on the assumptions of the model calculations), 
suffer the TDU that increases the C abundance in the surface trans¬ 
forming the O-rich star into a C-rich star; these stars have C/O > 1. 
In the most massive progenitors, with masses above ~ 4 - 5 Mg, 
HBB counteracts the effect of the TDU and prevents the formation 
of a C-rich star; therefore, for these stars one gets again C/O < 1. 

Moreover, the value of C/O in the atmospheres of AGB stars, 
which are amo ng the most efficient sources of dust in the Galaxy 
(IWhittelll201^ . defines the type of dust grains that are formed. 
When oxygen is more abundant than carbon, oxygen-rich grains, 
such as silicates and oxides, will be formed. If carbon is more abun¬ 


dant, other dust species are expected to form, such as SiC, TiC, and 
MgS. If we assume that the C/O abundance ratio (where C and O 
are the total abundances of these elements present in dust grains 
and in the gas) has not changed from the time when dust forma¬ 
tion initiated, the type of grains found in PNe will be related to the 
value of C/O in the atmosphere of the PN progenitor. Therefore, 
using infrared spectra we can identify different dust grains in PNe 
and thus classify them as C-rich (C/O > 1) or O-rich (C/O < 1) 
PNe. This inference on whether the ejected atmosphere of the pro¬ 
genitor star has C/O above or below 1 can be better than the one ob¬ 
tained with the C/O abundance ratio derived using emission lines, 
since the latter can be affected by depletion of C or O into dust 
grains and by the uncertainties invol ved in its calculation (see e.g. 
IPelgado-Inglada & Rodn'gue3l2014l) . 

Our aim here is to study the chemical abundances of a group 
of PNe that have high quality optical and infrared spectra. The op¬ 
tical spectra are used to derive the element abundances and the in¬ 
frared spectra are used to identify PNe with either oxygen-rich dust 
(ORD) or carbon-rich dust (CRD). We explore if the two groups 
of PNe, characterized by different dust grains, show differences in 
their chemical content. The abundances of the PNe are also com¬ 
pared with the predictions of recent nucleosynthesis models and 
with the chemical abundances derived for a group of Galactic H II 
regions. 


2 THE SAMPLE 


For the purpose of the present study, we need a sample of PNe with 
high quality spectra (i.e., with a high signal-to-noise ratio and spec¬ 
tral resolution better than 4 A), allowing us to perform the usual 
plasma diagnostics with good precision, to derive accurate abun¬ 
dances of elements that show only weak lines, such as Cl, and to 
determine (in most cases) the carbon-to-oxygen abundance ratio 
from very weak r ecombination lines. Most of the 57 pla netary neb¬ 
ulae compiled bv IPelgado-Inglada & R^rfg uez ( 20141) meet these 
requir ements. From this sample. IPelgado-Inglada & RodrfguezI 
1 I 2 OI 4 I 1 extracted 33 PNe with available infrared spectra. In this pa¬ 
per, we have selected the 20 PNe that we can classify as oxygen- 
rich or carbon-rich, i.e., objects that have in their infrared spec¬ 
tra either ORD features (amorphous or crystalline silicates) or the 
broad features around 11 and 30 yum associated with SiC and MgS, 
respectively, th at are expected to arise in carbon-rich environments 
(IWhittelll2003h . Many of the objects in the sample also show fea¬ 
tures from polycyclic aromatic hydrocarbons (PAHs), but since 
PAHs are sometimes detected in PNe with silicates, we do not 
use these features to classify the PNe as carbon-rich or oxygen- 
rich. Our selection criterion excludes NGC 2392, that contains dust 
grains (as can be inferred from the shape of its continuum), but 
shows no clear evidence of any of the previously mentioned dust 
features. We also excluded from the analysis PNe like NGC 6884 
and NGC 6741 that show PAH emission, but whose available in¬ 
frared spectra do not allow us to detect or to rule out the presence 
of silicates or the 11 and 30 yum features. 

In summary, we use the above information to classify the 
sample in two groups: 1) ORD PNe are those with silicates 
and 2) CRD PNe are those with the broad features associated 
with SiC and MgS. Some of the PNe with ORD and CRD also 
show emission from PAHs. The so-called dual chemistry PNe are 
those having PAH emission and silicate dust features (see e.g. 
IWaters et al.lll998l) . The compilation of dust features is taken from 
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Table 1. Sample of planetary nebulae 


Object 

Aliph. C 

Arom. C 

Silic. 

Comments 

Ref. 


0-rich dust PNe 

Cn 1-5 

No 

Yes 

Yes 

[WR] 

(1,7) 

DdDm 1 

No 

No 

Yes 


(2) 


H 1-50 

No 

Yes 

Yes 


(1) 


M 1-42 

No 

Yes 

Yes 


(3) 


M2-27 

No 

Yes 

Yes 

wels 

(1,8) 

M2-31 

No 

Yes 

Yes 

[WR] 

(1,7) 

M2-42 

No 

No 

Yes 

wels 

(1,9) 

MyCn 18 

No 

Yes 

Yes 


(4) 


NGC 3132 

No 

No 

Yes 

binary 

(4, 10) 

NGC 6210 

No 

No 

Yes 


(5) 


NGC 6439 

No 

Yes 

Yes 


(1) 


NGC 6543 

No 

No 

Yes 

wels 

(5,7) 

NGC 7026 

No 

Yes 

Yes 

[WR] 

(2, 7) 

C-rich dust PNe 

Hu 2-1 

Yes 

Yes 

No 


(2) 


IC418 

Yes 

Yes 

No 


(6) 


M 1-20 

Yes 

Yes 

No 

wels 

(1,11) 

NGC 40 

Yes 

Yes 

No 

[WR] 

(5,7) 

NGC 3242 

Yes 

No 

No 


(4) 


NGC 3918 

Yes 

No 

No 


(4) 


NGC 6826 

Yes 

No 

No 

binary 

(5, 12) 

References 

(DlWans & LiiH2007h. f21 

IWesson. Liu. & Barlow 

|2005 

), 

(3) 

Liu et a 

iJ J200lh. ( 4 ) IXsamis et al 

d2003h. (5) iLiu et alJ 

2004a 

). 

(6) 

Sharoee et al J2003b. 

17) iTvlenda. 

Acker. & Stenholml dl993b. 18) 

lOomv et al 

1 J2004h. ('9) iDeuew et a] 

J J2011 

). GO) ICiardullo et al.l 


[WR]: Wolf Rayet central star, wels: weak emission line 


IPelgado-Inelada & Rodrieuej ( 12014) and is based on Spitzer and 
ISO spectra. 

The sample of PNe includes objects from the Galactic disc, 
bulge, a nd halo. DdDm 1, with a height above th e Galactic plane of 
~ 7 kpc dOuireza, Rocha-Pinto. & Maciel|[200% . is our only object 
from the halo. Eight PNe belong to the bulge; Cn 1-5, H 1-50, M 1- 
42, M 2-27, M 2-31, M 2-42, N GC 6439, and M 1-20 These PNe 
satisfy the three criteria used bv IStasinska et al.l l ll99lh to identify 
bulge PNe, namely, they are located within 10° in terms of latitude 
and longitude with respect to the Galactic centre, their angular di¬ 
ameters are smaller than 20", and their radio fluxes at 6 cm are 
lower than 100 mjy. The remaining 11 PNe bel ong to the disc , and 
a calculation using the distances provided by IZhan 3 ( Il995h and 
adopting a solar Galactocentric distance of 8.0 kpc indicates that 
they have Galactocentric distances between 6 and 9 kpc, so that 
they can be loosely described as belonging to the solar neighbour¬ 
hood. 

Table [T] lists the 20 PNe studied here with the dust features 
identified in each of them, some comments about their central stars, 
and the references from where the optical line fluxes and the char¬ 
acteristics of the central stars have been taken. 

In addition to these PNe, we selected from the literature 
the seven Galactic H II regions from the solar neighbour¬ 
hood with the best available optical spect ra: M8, Ml6, M17 , 
M20, M42, NGC 3576 and NGC 3603 dEsteban et al.l |2004 
iGarcIa-Roias et alj2004Ll^006Ll2007h . They cover a range of Galac¬ 


tocentric distances that goes from 6.3 kpc (M16) to 8.65 kpc 
(NGC 3603), similar to the range covered by the disc PNe. 


3 THE ANALYSIS 
3.1 Physical conditions 

We use the software PyNeb l lLuridiana. Morisset. & Shawl2014h to 
calculate the physical conditions, the ionic abundances from col- 
lisionally excited lines (CELs), and the ionic abundances of He^ 
and He^^ from recombination lines (RLs). We have selected for 
our calculations the atomic data listed in Table 

We determine two electron temperatures, T^, for each nebula 
using the diagnostic ratios [O III] /l4363/(/l4959-l-/15007) and [N II] 
/15755/(/16548-l-/16583). We also compute an average density from 
the available ratios among the following: [O II] /13726//13729, [S II] 
T6716/T6731, [Cl ill] T5518/T5538, and [Ar IV] T4711/T4740. For 
each nebula, the adopted density is the median of the distribution of 
average values obtained in the Monte Carlo calculations that we use 
to estimate the uncertainties (see Section [T4t . We prefer to use the 
median since, in the low and high density limits, it is less affected 
than the mean by the large changes in the density caused by small 
variations in the intensity ratios. 

When is an important ionization state in a nebula, the re¬ 
combination of this ion may contribu te significantly to the intensity 
of the [N II] T5755 line (lRubinll986h . If this effect is not taken into 
accoun t, the value of 7’e([N II]) can be overestimated. iLiu et al.l 
( l2000h derived an expression to calculate this contribution that de¬ 
pends on the N^^ abundance, and also, to a small degree, on the 
value of Te- The N++ abundance can be determined using optical 
RLs, ultraviolet CELs, or infrared CE Ls, with the cho ice of lines 
leading to very different results (see e.g. lLiu et al .120001) . Therefore, 
the size of the correction is somewhat uncertain, and we decided not 
to use it here. However, we estimate an upper limit to this effect by 
using the largest possible value for the N^^ abundance, the one im¬ 
plied by RLs. The most affected objects are the PNe DdDm 1, NGC 
3242, and NGC 6826. In the case of DdDm I, where lWesson et ahl 
) l2005h estimate that recombination is the dominant contribution to 
the emission of [N II] 45755, we could use 7’e([0 III]) to derive all 
the chemical abundances, but we find very similar results. For NGC 
3242 and NGC 6826, the corrected values of TeljN II]) increase 
the O/H and Cl/O abundance ratios by ~ 0.2-0.6 dex. For the other 
objects, the corrections are much smaller, leading to changes in our 
abundance ratios lower than 0.1 dex. These changes, even the larger 
ones, do not affect in a significant way the conclusions we get from 
our results below and we will not consider them in what follows. 

Table [3] presents the derived temperatures and densities that 
will be used to compute the ionic abundances. 


3.2 Ionic abundances 

The ionic abundances of Ar^^, CE^, N^, Ne^^, O^, and are 
computed using the following CELs: [Ar III] 447136,7751, [Cl III] 
445518, 5538, [N ii] 446548, 6584, [Ne ill] 443869, 3968, [O ii] 
443726, 3729, and [O III] 444959, 5007, respectively. We do not 
consider the sulfur ions in our abundance determination because 
of the still unsolved “sulfur problem”, whereby the values derived 
in many PNe are found to be syste matically lower th an those in 
H II regions of the same metallicity dHenry et al .11 20041) . We adopt 
re([N II]) for the calculations involving the single ionized ions and 
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Table 2. Atomic Data 


Ion 

Transition Probabilities 

Collisional Strengths 

0+ 

0++ 

Zeinnen (1982) 

Storev & Zeinnen (2000i 

Kisielius et al. 120091 

Storev. Sochi. & Badnell (2014) 

N+ 

Galavis. Mendoza. & Zeinnen ('1997) 

TaYal(201D 

s+ 

S++ 

Mendoza & Zeinnen ('1982) 

Podovedova. Kelleher. & Wiese (2009) 

Taval & Zatsarinnv ("2010) 

Taval & Gupta (1999) 

Ne++ 

Galavis et al. Q997I 

McLauehlin & Bell 120001 

Ar++ 

Ar+2 

Munoz Buraos et al. 120091 

Mendoza & Zeippen (1982) 

Munoz Buraos et al. 120091 
Ramsbottom. Bell. & Keenan (1997) 

C1++ 

Mendoza Sz Zeinnen ('19!^2) 

Knieaer & Czvzak 119701 


Table 3. Physical conditions 


Object «e([0 It]) ne([S II]) «e([Cl III]) «e([Ar IV]) ne(adopted) TeCiN II]) TeClO III]) 
(cm“^) (cm“^) (cm“^) (cm“^) (cm“^) (K) (K) 


O-rich dust PNe 

Cn 1-5 . 

DdDm 1 4800+1“° 

H 1-50 6700+2“° 

M 1-42 1360+2“ 

M 2-27 . 

M2-31 . 

M 2-42 . 

MyCn 18 . 

NGC3132 . 

NGC 6210 . 

NGC 6439 3740+^“° 

NGC 6543 4500+1“° 

NGC 7026 2910+22 


^^^^-1100 

3000!2200 

7500'''^™ 

/ JOU_24oo 

8900+“°° 

izju_3oo 

1520+^‘^ 

7000'^^^®^ 

/ooo_2,00 

12600^^^®^ 

izouu_4g^^ 

6100+2500 

-1500 

O2?uU-2500 

3420+ 

■^^^^-850 

32OO+26O0 

•^^^Wi700 

-1300 

9200^^^^^ 

5?zuu_j200 

590+170 

-180 

740+350 

-280 

4070^^"^^^ 

3240+12°° 

4700+^^^^ 

^'^^-1100 

5060+q 

5900+^200 

-2100 

4400!«oo 

26500^^^®^^ 

8800!!“° 


'79Qn+2300 

2280 

3600^ 

J‘+^^_1200 

jyvv_24QQ 

‘tuuujgoo 

-1500 

9300“^^^^^ 

2?JUU_jgoo 

570+520 

j/U-310 

1240±220 

10700^^^®^ 

iu/uu_2400 

4190+900 

-710 


4100+4000 

HiuU-2goo 

-1200 

7400“^ 

-1100 

540+420 

-’‘+^-360 

621+200 
—o3 


38101™ 

6070+130° 

4980!°™ 

3400!|0° 

5100!™° 

5210+1100 

-710 

ouzu 7go 


7400+1^0 

/‘+«JU_2oo 

8680![“ 


12070±320 

11280!^“ 

10870!™ 

8860!™ 

9140±160 

9240!2™ 

8180!|“ 

11380!™ 

9760!|“ 

10060!™ 

8410!13° 

10020!™ 

7310+100 

'■^^^-120 

9350+280 

^■^^^-220 

9450+200 

11750+200 
/ ^^-400 

9970±290 

9630‘^200 

10260!™ 

9930+4^^ 

7780!™ 

9730+200 

-450 

9130!13° 


C-rich dust PNe 


Hu 2-1 

!:#UUU-2600 

IC418 

9000 

yuuu_23oo 

M 1-20 

^^^^-2700 

NGC 40 

1200+200 

NGC 3242 


NGC 3918 


NGC 6826 

1720+220 
^ ^^^-310 


28000!™°° 

15300!‘“°° 

10400!!™° 

9300!™° 

87oo!i?r 

1840!™ 

890±400 

2350'^^^^ 

zjjU-520 

1080±420 

5000'^^^®^ 

JUUU_J3QQ 

5480!1!«° 

1980+240 

1310™° 



11700!™° 

4860!g° 

9800™°° 

10600!“°° 

10900+"^^^^ 

L\jy\jyj_2gQQ 


1320!™ 

2100+”^^^ 

^^^+'-670 

I860!™ 

5920+^200 

5570!°™ 

2140+210 

1830!™ 


11790+240 

ii/vu_g7o 

9600!™ 

9530+200 

yjju_370 

8780!|“ 

11110!™ 

9760+2^^ 

^/UW-igo 

8460!^° 

10390±160 

12170+420 

'^-410 

11710+280 

il/lU-230 

10830!™ 

12540+200 

10400!“° 

9200!|“ 


H II regions 


M8 

1420+560 

1610+220 

1600+2 
iuuu_2gQ 

3900+2100 
”^^^^-2900 

1830!™° 

8290+‘40 

8040!™° 

M16 

1010+180 

lUlW-igo 

1420+200 

i‘+zU-27o 

1090+460 

i+'^^_490 


1180!™ 

8310!™ 

7580±130 

M17 

450+110 

520+160 

'^^^-140 

330+240 

djW-220 

2100!™° 

490!?!° 

8800!™ 

7960+ 

/!^dU-j3o 

M20 

230±50 

350±100 

390!™ 


310!*™ 

8250! 1“ 

7760!°° 

M42 

4930+1^00 

^^■^^-840 

6100+4200 

t>i^Wi800 

5790!™ 

4860!°°° 

5600+‘200 

jouU-680 

10200!™ 

8290!|™ 

NGC 3576 

1470+220 

13 30+200 

iJDU_32o 

2430!™° 

3000+^600 

duwu_j3oo 

2090+220 

8650!?™ 

8430±40 

NGC 3603 

2260+420 

4^o^480 

4000+^^®^ 

HUWU-iooo 

3800!!“° 

1700+2 
' -1100 

3160+660 

-510 

11120+460 

9010+130 


re([0 III]) for the double ionized ions; the adopted value of is 
used in all cases. 


The Ar++ abundance of DdDm 1 i s calculated using the 
[Ar III] /1/17135, 7751 line intensities of Henry et ^ ( 1200^ in¬ 
stead of those measured bv lWesson et al. I ll2005l) . the reference that 
we use for all the other line intens ities for this object. T he reason 
is that the intensities provided by IWesson et al. I J2005h for these 
lines are too low whe n comp a red to measurements by other authors 
([Barker & Cudwor ti] 1 19841 ; IClegg. Peimbert. & Torres-Peimberll 
1 19871 ; Henry et alT 2008h . which is unexpected since the angu¬ 


lar diameter of this object (0.6”) is smaller than t he slit widths 
used f or the observations. The intensities given by [Wesson et al. I 
( [2005h imply 12 -H log(A r++/H+) = 4.99, whereas the values listed 
by B^ker & Cudwortlj ( [l984l) ; [Clegg et ^ ( [l987[) ; [Henry et ^ 
( [2008[) cover the range 5.50-5.69. 


We use re([0 III]) and the average to compute the abun¬ 
dances of He+, He++, C++, and 0++from RLs. We calculate He+ as 
the average of the abundances derived from the He I lines 44471, 
45876, and 46678. The only exception is DdDm 1, where we only 
use the hrst two lines because the abundance of He+ derived from 
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A661S is too low. Since the intensity reported by IWesson et al. I 
( l2005h for this line is als o too low when compared to the mea¬ 
surements of other autho rs dSarker & Cudwo^ll984l : lciegg et al. I 
ll987l:lHenrv' et a lj2008ll. we think that there is some problem with 
the intensities of Wesson et al. I (l2005h for several lines in the red 
part of the spectrum, such as [Ar III] /1/17135,7751, He I /16678, 
and [S II] /16716. The He^^ abundances are derived using the 
H e II /14686 line. We use th e effective recombination coefficients 
ofjStore^^^ummra (|l^9^ for H I and He II and those computed 
bv IPorter et alj <2012112^ for He I. 

The abundances of are computed using C II 

/14267 and the case B effectiv e recombination coefficients of 
IPavev. Storey, & KisieliusI 1 I2OOOI) . As for the O^^abundances, we 
use the total intensity of multiplet 1 of O II, after correcting 
for the contribution of un d etecte d lines with the formulae of 
IPeimbert. Peimbert. & Rui3 (l2005h . and the recombination coef- 
ficients of istorevi ( Il993) . We exclude from our calculations O II 
lines of multiplet 1 reported as blended in the papers from which 
we took the line intensities. In the case of NGC 40, there is only 
one detected O II feature from this multiplet and it is blended with 
a N III line, thus w e prefer to use the abundance of derived 
bv iLiu et ^ (l2004bh from other multiplets. For the rest of the PNe, 
we do not consider other multiplets because their lines are weaker 
and likely suffer from blends. Indeed, we find that the abun¬ 
dances computed only from multiplet 1 are generally lower than 
those obtained using other multiplets. 

The final ionic abundances derived from CELs and RLs and 
their associated uncertainties (see Section [T4l are presented in Ta- 
bleH 


3.3 Total abundances 

To obtain total elemental abundances one needs to take into account 
the contribution of unobserved ions by using ICFs. Many studies 
still use ICFs based on similarities of ionization potentials. Those 
however ignore the fact that a small difference in the ionization po- 
tenti al may lead to a significant change in the ionization structur e 

(e.e. lSimon-DIaz & Stasihskal2008l : lDelgado-Inglada et ^l2Q14h . 

and that the ionization structure is not only governed by the spec¬ 
tral energy distribution of the stellar radiation field but also by 
the importance of recombination and charge exchange reactions. 
Photoionization models allow a much better estimate of ICFs. 
For PNe, the ICFs proposed bv iKingsburgh & Barlow! ( Il994ll have 
been used for two decades. However, there is no real documenta¬ 
tion on how they were obtained and there is no estimate of their 
uncertainties. A set of formu lae has been proposed recently by 
IPelgado-Inglada et al. I ( l2014h . based on an extensive grid of ab ini¬ 
tio photoionization models aiming at representing the whole man¬ 
ifold of PNe. The large set of considered photoionization models 
allowed, for the first time, an estimation of the error bars on ICFs. 
They are computed us i ng the analytical expressions provided by 
IPelgado-Inglada et al. I (|2014|) . which are based on the maximum 
dispersion of each ICF obtained from the grid of photoionization 
models. Those ICFs were tested by applying them to a large sam¬ 
ple of PNe and checking that the derived abundance ratios did not 
depend on the degree of ionization of the objects. 

The sample of PNe considered here is a specific sample ex¬ 
cluding high density objects and objects io nized by very hot stars. 
We th erefore checked if the ICFs from IPelgado-Inglada et al. I 
(|2014|) . which were obtained by considering a wider sample of 
physical conditions, remain valid for the restricted subsample of 
photoionization models better representing our observational sam¬ 


ple, i.e., we discard the models with very high gas density or with 
ionizing stars of very high effective temperatures. We found that 
the ICFs presented in Delgado-Inglada (2014) remain valid for the 
restricted subsample of models, both as regards the analytical for¬ 
mulae for the ICFs and as regards the error estimates. _ 

W e thus adopted the ICFs derived bv IPelgado-Inglada et al. I 
( I2OI4I) to calculate the abundance ratios Ar/0, C/0, Cl/0, Ne/O, 
He/H, and 0/H. The ICF for O is based on He++/(He+-l-He++), the 
ICFs for Ar, C, Cl, and He are based on 0++/(0^-l-0+^), and the 
ICF for Ne is based on both ratios. As fo r N, we show here the 
results obtained with the ICF derived by IPelgado-Inglada et al. I 
( I2OI4I) as well as those derived with the widely used ICF, N/O = 
N^/0^ (which are the ones listed in Table |3. We will discuss the 
differences between both estimates of N/O in section 4.5. 

For the H II regions we are considering in this paper, the situ¬ 
ation is a priori worse. Formula e for ICFs ba s ed on photoionization 
models have been proposed by lizotov et"^ (l200bh . but they were 
proposed and tested for giant extragalactic H II regions, and there 
is a priori no reason why they should apply to slit observations of 
nearby dense H II regions, some of them ionized by single stars. 
In fact, we find that they introduce a dependence of the total abun¬ 
dances of argon and neon on the degree of ionization for the H II 
regions in our sample. In absence of any adequate set of ICFs, we 
used for our H II regions the same expressions that we used for our 
PNe. In principle, we could trim the grid of models built for PNe 
to retain only those that would be adequate for the H II regions we 
consider (ie. remove the models with effective temperatures above 
50000 K and all those that have a shell structure or are density 
bounded) and use the trimmed grid to compute ICFs that would be 
relevant for our sampl e of H II regions. In practise, w e checked that 
the formulae given bv IPeleado-Inglada et al. I ( 1201 4l) actually fit the 
trimmed grid quite well (except for the determination of Ar/0 in 
objects where 0^^/(0^-l-0^^) is larger than 0.5, where the resulting 
bias on Ar/O can go up to ±0.09 dex. 

One important advantage of the ICFs used here to determine 
the Cl and Ar abundances is that they can be used when only lines 
of CP^ and Ar^^, respectively, are observed. These ions can be 
analysed using relatively bright lines in the optical spectrum, which 
makes it easier to achieve a homogeneous calculation of the total 
abundances. Other ICFs in the literature include more ions that have 
relatively weak emission lines, such as the [Ar IV] /I/I4711,4740 
lines in some objects, or that lie in a spectral range that it is not often 
observed, such as [Cl II] A9122. Moreov er, we found, from the 
:rid of photoionization models computed in lPelgado-Inglada et al. I 
20141) . that the other ICFs in the literature for Ar and Cl introduce 
a higher dispersion in the derived abundances. 

As for the 0/H values comp uted from RLs , for th e H II re¬ 
gions we adopt those provided by lEsteban et al.l ( l2005h . They are 
based on the and abundances derived with O I and O II RLs 
for some objects; for other objects O I RLs could not be measured 
properly and the abundances are based on the abundances im¬ 
plied by CELs corrected upwards to account for the possible pres¬ 
ence of small tem perature fluctuations within the nebulae. We note 
that the results of lEsteban et al.l ( l2005h show that some of the dis¬ 
persion in these values of O/H (and in the corresponding values 
based on CELs) arises from the Galactic abundance gradient. The 
assumpti on of small temperat ure fluctuations might not be right for 
the PNe dWesson et al. 11200^. Besides, only I C 418 has measure¬ 
ments of suitable O I RLs 1 Sharpee et alj2003h . Hence, we derived 
their RLs oxygen abundances using the abundances implied by 
multiplet 1 of O II and the ionization fraction for this ion inferred 
from the CELs analysis. 












































































6 Delgado-Inglada et al. 


Table 4. Ionic abundances from CELs and RLs: X*' = 12 + log(X+7H^) 


Object 

Ar++ 

C1++ 

N+ 

Ne++ 

0 + 

0 ++ 

He+ 

He++ 

C++ 

0 ++ 


CELs 

CELs 

CELs 

CELs 

CELs 

CELs 

RLs 

RLs 

RLs 

RLs 

O-rich dust PNe 

Cn 1-5 

6.38±0.03 

5.22±0.05 

8 19+0.07 
°-^^-0.03 

8.28±0.04 


8.70±0.04 

11 . 01 ± 0.01 


9.09 

9.01 

DdDm 1 

5.61±0.04 

4.59±0.04 

6 75+0.08 
^•^'^-0.04 

7.19±0.06 

7 'JS+0.20 

'•■^°-0.09 

7.87±0.04 

10.94±0.02 



8.39 

H 1-50 

5.99±0.03 

4.88±0.04 

7 01+0.06 
-0.03 

7.94±0.04 

7 'JQ+0.13 

8.64±0.05 

11 . 00 ± 0.01 

10.04±0.02 

8.35 

9.09 

M 1-42 

6.27±0.03 

5.04±0.04 

7.89±0.04 

8.03±0.04 

-j 64+0.07 
'•°^-0.04 

8.39±0.04 

11 . 20 ± 0.01 

10.04±0.02 

9.39 

9.63 

M2-27 

6.50±0.03 

c 9Q+0.06 
'^■^^- 0.02 

7.66±0.07 

8.37±0.04 


8.84±0.04 

11.15±0.01 

8.84±0.02 

8.85 

9.27 

M2-31 

6.15±0.04 

5.00±0.04 

7.04±0.05 

8.06±0.04 

'•^^-0.06 

8.64±0.05 

11.07 ±0.01 




M2-42 

6.16±0.03 

5.14±0.07 

' -^^-0.03 

8.04±0.04 

7 "ja+O.n 
'•-’-’-0.06 

8.73±0.04 

11.06±0.01 

8.45±0.08 



MyCn 18 

6.33±0.03 

5.31 ±0.03 

7.57±0.04 

8.18±0.05 

7 61 S+O.O 9 
'■°°-0.06 

8.53±0.04 

11 . 00 ± 0.02 

8 . 66 ± 0.02 

8.36 

8.87 

NGC3132 

6.38±0.03 

5.18±0.03 

8.19±0.04 

8.30±0.04 

8.47±0.06 

8.53±0.05 

11.05±0.01 

9.52±0.02 

8.83 

8.81 

NGC 6210 

5.93±0.04 

4.70±0.05 

6.47±0.05 

8.12±0.06 

7.09±0.10 

8.55±0.06 

11 . 02 ± 0.02 

9.28±0.02 

8.80 

8.97 

NGC 6439 

6.35±0.03 

5.12±0.03 

7.52±0.05 

8.16±0.04 

7.67±0.07 

8.60±0.04 

11.05±0.02 

10.31±0.02 

8.99 

9.06 

NGC 6543 

6.48±0.03 

5.22±0.06 

6.80ti? 

8.24±0.05 

'•^^-0.06 

8.79±0.05 

11.06±0.01 


8.77 

9.07 

NGC 7026 

6 . 22 ± 0.02 

5 16+0.06 


8.19±0.03 

'•^^-0.06 

8.64±0.03 

11.03±0.01 

10 . 11 ± 0.02 

8.94 

9.07 

C-rich dust PNe 

Hu 2-1 

5.69±0.03 

A C 1 +0.07 
-0.03 

6 93+0.16 
^•^'^-0.04 

7.44±0.04 

7 

'•^^- 0.11 

Q 7 - 3 + 0.07 
°-"^^-0.04 

10.91±0.01 

8.39±0.02 

8.62 

8.71 

IC418 

6.00±0.03 

4.82±0.04 

7.63±0.06 

6.78±0.04 

°-^^-0,15 

8.09±0.04 

10.96±0.01 


8.74 

8.24 

M 1-20 

5.77±0.03 

^•^^-0.06 

6 79+0.11 
^•'^-0.05 

7.69±0.04 

7 34 + 0.28 
'•■^^- 0.10 

8.55±0.05 

10.99±0.01 

'•^^-0.16 

8.67 

8.66 

NGC 40 

5.67±0.02 

4.65±0.03 

7.94±0.04 

5.80±0.03 

8.65±0.06 

7.11 ±0.03 

10.78±0.01 

-0.16 

8.81 

8.33 

NGC 3242 

5.60±0.03 

4.40±0.03 

5.55±0.04 

7.88±0.04 

6.50±0.06 

8.43±0.05 

10.87±0.02 

10.34±0.02 

8.79 

8.84 

NGC 3918 

5.93±0.03 

4.71 ±0.04 

7.15±0.05 

7.88±0.04 

'•^^-0,06 

8.45±0.05 

10.82±0.02 

10.55±0.02 

8.70 

8.82 

NGC 6826 

6.06±0.02 

4.85±0.03 

6.22±0.06 

7.88±0.03 

7.00±0.08 

8.55±0.04 

11 . 00 ± 0.01 

7 34 + 0.10 
-0.16 

8.74 

8.82 

H II regions 

M 8 

6.16±0.02 

5.04±0.03 

7.54±0.03 

7.05±0.04 

Q - 5 . 7 + 0.12 
-0.04 

7.89±0.03 

10.84±0.01 


8.31 

8.24 

M16 

6.19±0.04 

5.09±0.04 

7.73±0.03 

7.15±0.05 

8.47±0.05 

7.92±0.04 

10.90±0.01 


8.39 

8.31 

M17 

6.30±0.03 

5.06±0.03 

6.84±0.04 

7.76±0.03 

7 80+0.07 
"■°^-0.04 

8.45±0.03 

10.98±0.01 


8.73 

8.72 

M20 

6.16±0.02 

5.04±0.03 

7.59±0.04 

6.68±0.03 

8.48±0.05 

7.73±0.02 

10.85±0.01 


8.19 

8.18 

M42 

6.33±0.03 

5.12±0.04 

6.87±0.04 

7.77±0.03 

7 71 + 0.10 
-0.05 

8.43±0.05 

10.98±0.01 


8.35 

8.61 

NGC 3576 

6.30±0.05 

4.94±0.03 

7.06±0.04 

7.70±0.02 

8.11 ±0.06 

8.37±0.01 

10.97±0.01 


8.44 

8.67 

NGC 3603 

6.29±0.02 

5.04±0.04 

6.48±0.06 

7.78±0.03 

7 34 + 0.12 
'•-’‘+-0.06 

8.44±0.03 

11 . 00 ± 0.01 


8.49 

8.71 


The final abundances for each PN and H II region are listed in 
Table |5] 


3.4 Uncertainties 

The uncertainties in the physical conditions and the ionic abun¬ 
dances are computed through a Monte Carlo simulation. We gen¬ 
erate 200 random values for each line intensity using a Gaussian 
distribution centred in the observed line intensity and with a sigma 
equal to the associated uncertainty. For higher number of Monte 
Carlo simulations, the errors in the computed quantities remain the 
same. 

We derive the values of T^, n^, and the ionic abundances for 
every Monte Carlo run and calculate the errors associated with each 
quantity from the 16 and 84 percentiles; these values define a con¬ 
fidence interval of 68 per cent which is equivalent to one standard 
deviation in a Gaussian distribution. We assign a typical uncertainty 
of 5% and 15% to the C^^ and ionic abundances, respectively, 
computed from RLs. We use this simplified approach because the 
determination of uncertainties for the abundances implied by 
RLs is difficult and unreliable, and the C^^ abundances are only 
used to calculate the C/O abundance ratio from the ratio C^^/O^’^. 


Some of the O II lines from multiple! 1 are not observed and we 
have to estimate their contribution to the multiple!, which can de¬ 
pend on the electron density. In addition, some of the O II lines 
could be blended with nearby weaker lines. An uncertainty of 15% 
corresponds to the typical difference between the abundances de¬ 
rived from individual lines of multiple! 1 and from the total in¬ 
tensity of the multiple! and al so to the differences with the 
compu ted with other multiplets l lGarcIa-Roias et alj2013l;lLiuetal.l 
l2004bh . A complete analysis of the uncertainties related to the 
abundances implied by RLs should also take into account the infor¬ 
mation provided by the other multiplets and would be rather oner¬ 
ous. We consider that our simple approach is sufficient for the pur¬ 
poses of this paper. The uncertainty of 5% for the abundances 
arises from the typical uncertainty of the intensity of C II 4267 
relative to HjS in our sample objects. 

As for the total abundances, we consider not only the uncer¬ 
tainties associated with the physical conditions and the ionic abun¬ 
dances but also those related to the adopted IGF. To estimate the 
uncertainties associated with the ICFs we perform another Monte 
Carlo simulation. We construct a uniform distribution for each IGF, 
where the central value is the IGF calculated from the observed de¬ 
gree of ionization, and the minimum and maximum values are ob- 
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Table S. Abundance ratios: {X/H| = 12 + log(X/H) 


Object 

(Ar/H) 

(Cl/H) 

(N/H)* 

(Ne/H) 

(O/H) 

I He/H) 

log(C/0) 

{O/H} 

Galactic 


CELs 

CELs 

CELs 

CELs 

CELs 

RLs 

RLs 

RLs 

component 

O-rich dust PNe 

Cn 1-5 

^ 4ff+0.13 

5.35±0.08 

0 77 + 0.12 
-0.24 

0 ^^+0.09 
®-™-0.04 

8.84!®:“ 

11 . 10 ± 0.01 

0 

-0.38 

9.14 

bulge 

DdDm 1 

5 69+0.14 
^•^^-0.32 

4 7T+0.04 
^•'•^- 0.10 

-0.25 

7.55±0.07 

7.99!®;®’ 

10.94±0.02 


8.51 

halo 

H 1-50 

°-^^-0.24 

c 1 0 +O.I 8 
^•^°- 0.12 

0 07 + 0.11 
-0.27 

7 96+0.03 
-0.07 

8.69±0.04 

11.04±0.01 

-0.68±0.38 

9.14 

bulge 

M 1-42 

6 42+0.15 
'^•^^-0.28 

c 7 T+0.03 
'^•"^^- 0.10 

8.73t>5 

8.05±0.07 

8.48±0.04 

11 . 22 ± 0.01 

^•^^-0.40 

9.72 

bulge 

M2-27 

^•^^-0.23 

-0.09 

8.89/«:>^ 

8.54±0.06 

8 Q7+0.06 

-0.03 

11.16±0.01 

-0.37±0.38 

9.30 

bulge 

M2-31 

-0.08 

c 01 +0.26 

-^■^^-O.IO 

8.46!«:>5 

8.20±0.06 

8.66±0.04 

11.07±0.01 



bulge 

M2-42 

6 37+0.56 
-0.16 

5 45+0.24 
■^■^'^- 0.12 


8.18±0.05 

8.74±0.04 

11.06±0.01 


9.50 

bulge 

MyCn 18 

^ 47+0.15 
-0.32 

- 0.11 

8.48!«:>^ 

8.41±0.06 

8.58±0.04 

11 . 00 ± 0.01 

-0.49±0.37 

8.93 

disc 

NGC3132 

ft, zin+0.14 
°-^^-0.27 

■^■^^- 0.10 

Q co+0.14 
O.JJ -020 

8.36±0.06 

8.81±0.03 

11.06±0.01 

-0.13±0.34 

9.08 

disc 

NGC 6210 

^•^'^- 0.10 

5 o 2 +^- 2 ^ 
'^■^^-0.06 

7 QC+0.15 
^•^^-0.24 

8.16±0.07 

8.57±0.05 

11 . 02 ± 0.02 

-0.33 

8.99 

disc 

NGC 6439 

^•^'^-0.25 

c ' 17 + 0.04 
- 0.10 

8.54!®;>5 

8.20±0.05 

8.69±0.05 

11 . 12 ± 0.02 

-0.05± 

9.15 

bulge 

NGC 6543 

6 70+0-60 
°-'^-0.08 

c C 7 +O .26 

J.P /-0 11 

8.40!®;^ 

8.36±0.07 

Q q/a+0.06 

-0.03 

11.06±0.02 

A 90 + 0.40 

^•^■^-0.36 

9.08 

disc 

NGC 7026 

°-^^-0.34 

c 41 +0.05 
-0.09 

8.62!®:>3 

8.22±0.05 

071 +0.06 
' - 0.02 

11.08±0.01 

-0.38 

9.14 

disc 

C-rich dust PNe 

Hu 2-1 

C Q 1 +0.14 
-0.31 

A ^^8+0.05 
^•°^-o.io 

7.73±0.16 

771 + 0.20 

-0.04 

8.30!°;^^ 

10.91±0.01 

-0 

^•^^-0.39 

8.78 

disc 

IC418 

^•^'^-0.32 

4 94+0.04 
- 0.11 

7.82/®:|® 

7 54+0.05 
- 0.12 

8.54±0.09 

10.96±0.01 

^•^'^-0.34 

8.69 

disc 

M 1-20 

5 96+0.21 

4 99+0.11 
- 0.12 

8-03tl6 

7 oc+ 0.11 

'•°^-0.05 

8 <q+0.06 

®-^°-0.03 

10.99±0.01 

0.06/®:®® 

8.68 

bulge 

NGC 40 

6 17+0-‘4 
-0.63 

4 99+0.05 
^•^^- 0.10 

7 96+0.15 

'■^^- 0.11 

-0.15 

8.66±0.05 

>10.78 

-0.48/®:®® 

9.89 

disc 

NGC 3242 

c Q 1+0.61 
-0.06 

4 09 + 0.25 
^•°^- 0.12 

- 0.12 

7 9 C+ 0.04 
'•^'^-0.08 

8.50±0.06 

10.98±0.01 

^•^■^-0.34 

8.91 

disc 

NGC 3918 

f. 1 O+0.15 
^•^°-0.36 

5 01 +0-05 

8.13±0.14 

8.00±0.08 

8.63±0.06 

11 . 01 ± 0.01 

-0.12±0.36 

9.01 

disc 

NGC 6826 

r 90 +O. 6 O 
^•^^-0.05 

c 19 + 0.26 
'^•*^-0.09 

7 7Q+0.16 

^■^°-o.io 

8.01±0.06 

8.56/®;®® 

11 . 00 ± 0.01 

_0 01+0-39 
-0.36 

8.83 

disc 

H II regions 

M 8 

6 

^•^“’-0.31 

c 1 O+0.06 

'^•10-0.09 

7 67+^-^^ 
- 0.11 

7 65+®-^^ 
'•^'^-0.05 

8 <n+o io 
-0.03 

>10.84 

-0.27±0.30 

8.71 

disc 

M16 

-0.25 

5 94+0.05 
-^■^^- 0.11 

7 g 4 + 0.18 
^•°^-o.ii 

7.81±0.05 

8.58±0.04 

>10.90 

-0.29±0.33 

8.80 

disc 

M17 

^ 41 + 0.15 
°-^^-0.34 

c 99+0.05 
'^■^^-0.09 

7.58±0.13 

7.85±0.05 

8.54±0.03 

10.98±0.01 

^•^^-0.33 

8.76 

disc 

M20 

6 32+^-^^ 

-0.31 

c 9 1 +0.04 
-^■^^-O.IO 

7.66!®;lf 

7.50±0.06 

8.55±0.04 

>10.85 

_n 47+0-34 
^•^'-0.29 

8.71 

disc 

M42 

6 46+^- 
°-^^-0.29 

c 99+0.04 
-^■^^-O.IO 

7-6811^ 

7.84±0.05 

8.51 ±0.05 

10.98±0.01 

-0.33 

8.65 

disc 

NGC 3576 

6 34+0.13 
-0.30 

■^■^^-0.10 

'■'^^-0.08 

7.89±0.05 

8.56±0.03 

10.97±0.01 

-0.33±0.34 

8.74 

disc 

NGC 3603 

6 47+0.14 
°-^^-0.30 

c 9 Q+U.UD 

D.28-0 11 

7.60±0.13 

7.82±0.05 

8.47±0.03 

11 . 00 ± 0.01 

-0.30 

8.72 

disc 


* N/H computed from N/0 = N+/0^. 


tained from the uncerta i nties a ssociated with each ICF provided in 
IPelgado-Inglada et al. I l l2014h . From this distribution we generate 
random values that we use to calculate the total abundances for the 
200 runs. The errors associated with the abundance of each element 
are given by the 16 and 84 percentiles. Our computed uncertainties 
in the total abundances are probably overestimated since we are us¬ 
ing the maximum dispersion in each ICF, created by a large grid of 
photoionization models, whereas real nebulae are probably better 
represented by a smaller set of models. 

In the case of nitrog en, we use the method describ ed above 
when using the ICF from IPelgado-Inglada et al. I l l20I4h . and we 
adopt a uniform distribution of errors of ±0.2 dex when using N/0 
= N^/O^, which seems reasonable given the dispersion we get for 
abundance ratios like N/0 and N/Cl in the sample H 11 regions. 

As for He, IPelgado-Inglada et al. I ( l2014h proposed that one 
should compute the total He abundance by adding the contribu¬ 
tions of He^ and He^^, including the contribution of He° in the 
uncertainties associated with He/H. There are four objects in the 
sample (three H II regions and one PN) with a very low degree of 
ionization, where the abundance of He® is important and thus, the 
value of He/H is very uncertain. We represent their He abundance 


in the figures as a lower limit and we only plot the errors associated 
with the He ionic abundances and n ot with the ICF. _ _ 

Each of the ICFs derived in IPelgado-Inglada et al. ] i2014h 
is valid for a certain range in the degree of ionization, given by 
0++/(0^+0^^) or He++/(He+±He++); outside this range the ICF 
formula is not necessarily c orrect. Here we use the ICFs from 
IPelgado-Inglada et al. I l l2014ll in all the objects, regardless of their 
degree of ionization, but we assign higher uncertainties in the ICFs 
when the objects fall outside the range of validity. The low disper¬ 
sion found in some of the abundance ratios, such as Ar/Cl, suggests 
that the adopted ICFs are still valid outside the ranges proposed in 
IPelgado-Inglada et al. I l l2014h . at least in some cases. 

Since the ICFs for many elements are based on the 
0 ^^/(0^+0^^) ratio, the uncertainties introduced by these correc¬ 
tions are not independent. If we do not consider this effect, we get 
very large uncertainties for some of the abundance ratios consid¬ 
ered be low, namely, Ar/Cl and Ne/Cl . We have used the model re¬ 
sults of IPelgado-Inglada et al. I l l2014ll to obtain better estimates of 
the uncertainties in those cases. The uncertainties in log(Ar/Cl) as¬ 
sociated with the ICF are: 

±0.30/ - 0.90 dex, for 0“"-"/(0+ ± 0++) < 0.05 
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+0.20/ - 0.15 dex, for 0.05 < 0++/(0+ + 0++) < 0.45 
+0.10/ -0.04 dex, for 0.45 < 0++/(0+ + 0++) < 0.95 
+0.66/ - 0.01 dex, for 0++/(0+ + 0++) > 0.95. 

As for the log(Ne/Cl) values, the uncertainties associated with the 
ICF are; 

+0.70/ - 1.30dex,forO""*/(O+ +0"^-") < 0.1 
+0.13/ -0.08 dex,for 0.1 < 0++/(0+ + O"'"’) < 0.95 
+0.20/ - 0.30 dex, for 0++/(0+ + > 0.95. 

We do not include the uncertainties on the atomic data since 
it is very difficult to quantify them hut we h ave selected the atomi c 
data that we consider more reliable (see e.g. lStasihska et alll2013h . 

The final adopted value for each parameter is the value com¬ 
puted from the observed intensities except for the electron density, 
for which we prefer the median of the distribution of average densi¬ 
ties obtained from the available diagnostics in each run. In the low- 
and high-density limits, the median is less affected by the large 
variations obtained in the density values due to small variations in 
the intensity ratios. Note that the electron density is a special case; 
for some of the other calculated quantities, the use of the median 
would lead to similar results, for other quantities, such as the total 
abundances, it is not a good idea to use the median since their dis¬ 
tributions depend on the distribution of the uncertainties due to the 
ICFs, for which we can only estimate the width, not the shape. 


4 RESULTS AND DISCUSSION 
4.1 Chlorine as a proxy for metallicity 

We have derived the abundances of He, C, N, O, Ne, Cl, and Ar. 
From these elements, only Cl and Ar are not expected to be modi¬ 
fied by the evolution of low and intermediate mass stars, and can be 
used to infer the metallicity of the PN progenitor star. We could also 
use S but, as we mentioned in Section 1, there are problems with the 
determination of the abundance of this element. Figure[T]shows the 
Ar/Cl abundance ratio as a function of Cl/H for the group of PNe 
and H II regions studied here. The small dots in the figure corre¬ 
spond to the values derived for each Monte Carlo run and take into 
account all the uncertainties arising from errors in the line fluxes in¬ 
volved in the calculations and those related to the adopted ICF. The 
metallicity of our objects changes by an order of magnitude, but 
all of them show log(Ar/Cl) - 1.1 with very low dispersion. The 
average value of the Ar/Cl abundance ratio considering our whole 
sample of ionized nebulae is log(Ar/Cl) = 1.10 ± 0.02, where 0.02 
is the standard error of the mean. The tight correlation between the 
Cl/H and the Ar/H values is in agreement with a lockstep evolution 
of both elements, at least for the metallicity range covered by our 
objects. Since the Ar and Cl ab undances are v ery difficult to calcu¬ 
late in the inte rstellar medium l|jenkinsll2009h . in stars, and in the 
Solar System jAsplund et alj|2009h . the results obtained from the 
analysis of photoionized nebulae provide the best estimates of the 
abundances of these two elements. 

The uncertainties in the values of Ar/H are significantly higher 
than those found for Cl/H, as can be seen in Table Hence, in 
what follows we use the Cl/H abundance ratio as a proxy for the 
metallicity of all the objects. 


2.0 



4.4 4.6 4.8 5.0 5.2 5.4 5.6 5.8 6.0 

12 + log(CI/H) 

Figure 1. Values of Ar/Cl as a function of Cl/H. The green circles represent 
ORD PNe, the red diamonds CRD PNe, and the blue squares the H II re¬ 
gions. We mai'k with open diamonds the PNe with PAH emission and with 
a small black dot the PNe from the bulge. The coloured small dots corre¬ 
spond to the values obtained for the 200 Monte Carlo run s performed for 
each object. The protosolar abundances of iLoddersI hoid) are overplotted 
with the solar symbol. The CRD PN NGC 6826, with 12 + log(Cl/H) = 5.19 
and log{Ar/Cl) = 1.09, is hidden behind the H II regions. 

4.2 Metallicity differences between H ii regions and PNe 
with CRD and ORD 

We also show i n Figure [T] the protosolar abundances provided by 
iLodderd ( l2010h . Unfortunately, the S olar System abundan ces of 
these elements are not well constrained dAsolund et alJ2009h . If we 
use the photospheric abundances provided bv lAsplund et al.l ( l2009t) 
for Ar and Cl, 12 + log(Ar/H) = 6.40 and 12 + log(Cl/H) = 5.50 
with a correction of ~ 0.04 dex due to diffusion effects, we would 
get 12 + log(Cl/H) = 5.54 a nd log (Ar/Cl) = 0.9 instead of the 
values favoured by lUoddersI ( l201oli . 12 + log(Cl/H) = 5.3 and 
log(Ar/Cl) = 1.25. However, it is clear from Figure Q] that most 
of the ORD PNe have abundances close to or above solar, H II re¬ 
gion abundances are below solar, and CRD PNe show even lower 
abundances (note that NGC 6826, our more metal-rich CRD PN, 
with 12 + log(Ar/H) = 6.28 and 12 + log(Cl/H) = 5.19, is hid¬ 
den behind the H II regions in this figure). It seems that the pro¬ 
genitors of the CRD PNe in our sample were formed in a sub¬ 
solar metallicity medium, which agrees both with the expectation 
of a higher efficiency of the TDU at low metallicities (see e.g. 
iKarakas & Lattanziol2014h and with the observation of larger frac¬ 
tions of CRD PNe in the metal-poo r environment of the Magellanic 
Clouds teernard-Salas et alj2009ti . On the other hand, the progeni¬ 
tors of the ORD PNe were born in env ironments that cover our full 
metall icity range. In agreement with iGarcfa-Hernandez & Gornvl 
( l2014h . we find that PNe with dual chemistry (with both silicates 
and PAHs in their spectra) have metallicities close to solar. 

Several ORD PNe (7 out of 13) belong to the bulge, where 
the aver age PN metallicity is known to be higher than in the disc 
(see e.g. lWang & LiuluOOTh . However, three of the disc PNe have 
higher values of Cl/H than the H II regions. Since all our disc ob¬ 
jects are located roughly at the same range of Galactocentric dis¬ 
tances, chemical evolution models predict that the PNe progenitor 
stars should have formed at lower metallicities than the ones traced 
by H II regions. iRodrfguez & Delgado-Ingladal l l201lh obtained a 
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Fi gure 2. Values o f Ne/Cl as a function of Cl/H. The protosolar abundances 
of lhodde^ j201(]h are overplotted with the solar symbol. We use the same 
symbols as in Figure[T] 


similar result from an analysis of the oxygen abundances in a group 
of PNe and H II regions belonging to the solar neighbourhood. 
iRodrfguez & Delgado-Ingladal l l201 ih speculated that this could be 
due to the presence of organic refractory dust grains in H II regions 
containing significant amounts of oxygen, but the fact that we find 
the same result for chlorine indicates that this is a more general ef¬ 
fect. A possible explanation is stellar migration of the central stars 
(and the Sun?) from inner regions of the Galaxy, but other expla¬ 
nations, such as infall of metal-poor gas clouds into the Galaxy or 
changes in the stellar composition arising during star formation or 
stellar evolution cannot be ruled out at the moment. Six of the seven 
H II regions show similar metallicities, indicating that the interstel¬ 
lar medium in the solar neighbourhood is chemically well mixed. 
This is confirmed by the values we find for O/H for the seven H II 
regions (see Tablej^and Figure[^below). The chlorine abundance 
of NGC 3576 is somewhat low, but compatible within the errors 
with the values of the other H II regions. 


4.3 Inconclusive evidence for neon enrichment 

Figurel^shows the results for neon. Neon ca n be produced in stars 
with i nitial masses in the r ange ~ 2-4 jKara kas & Lattanzid 
I 2 OO 3 I : IPignatari et alj|201^ . and iMilingo et alj l l2010t) argue that 
there is evidence for neon enrichment in some Galactic PNe. Our 
Ne/Cl abundance ratios show a large dispersion, and the fact that 
CRD PNe, FI II regions, and ORD PNe cover similar ranges in 
Ne/Cl suggests that this dispersion arises from uncertainties in the 
ICF for neon. Our sample PNe do not show clear evidence for Ne 
enrichment, although we cannot rule it out since it could be masked 
by the uncertainties in the derived Ne/H abundance ratio. If we re¬ 
ject the high value of Ne/Cl in NGC 3576, since it could be due to 
the too low value found for Cl/H in this object, we can see that the 
other H II regions and several PNe cluster around the solar value of 
Ne/Cl. The remaining PNe might have some neon enrichment. 



Figure 4. Values of 0/Cl as a function of 0^^/(0^+0^^). The protosolai' 
abundances of lLodder^ j2010li are overplotted with the solar symbol. We 
use the same symbols as in Figure[T] 


4.4 Oxygen production in the PNe with CRD 

Our results for the O/Cl abundance ratio as a function of O/H and 
Cl/H are presented in Figure Note that the trend followed by the 
dots that show the results of our Montecarlo propagation of errors 
in the right panel of Figure]^ is caused by the correlation between 
the uncertainties in the values we derive for 0/Cl and Cl/H. The 
trend is not present in our computed values for these abundance 
ratios. However, in order to check if the ICF used to calculate the 
O/Cl values is introducing a bias in our results, we plot in figure |4] 
the O/Cl abundance ratios as a function of the degree of ionization, 
given by m = 0^^/(0^ - 1 - O^^), which is the parameter on which the 
ICF is based: Cl/O = (CF+/O+)(4.1620 - 4.1622m'’'2‘)“’5. We see 
that the adopted ICF does not seem to be introducing any kind of 
bias in our results. In particular, the CRD PNe cover the full range 
of values of ai and their O/Cl abundance ratios do not follow any 
trend with ai. 

It can be seen in Figure that the metallicities implied by the 
Cl and O abundances do not agree for all the objects. The Sun, 
the H II regions, and the ORD PNe show similar values of the O/Cl 
abundance ratio, in agreement within the errors. However, the CRD 
PNe have values of O/Cl which are clearly higher. The exception 
is NGC 6826, the CRD PN with the highest metallicity, according 
to its chlorine abundance. The straightforward explanation of this 
result is that the low-metallicity CRD PNe besides being C-rich are 
enriched in oxygen. 

As we commented in Section 1, canonical AGB models like 
those of iKarak^ ( 1201 Oh do not predict oxygen enrichment at the 
metallicities that characterize our sample P Ne, but models with dif - 
fusive convective overshooting like those of IPignatari et al.l (1201 3l) . 
predict that the ejecta of stars with initial masses of ~ 2-3 Mq and 
metallicities Z = 0.01 and 0.02 will be enriched in both C and O 
(with C/O >1), with an increase in the O/H abundance ratio of up 
to ~ 0.26 dex. This agrees with the differences we find in the av¬ 
erage values of log(0/Cl) for our objects: log(0/Cl) = 3.65 ± 0.06 
for the CRD PNe excluding NGC 6826, 3.34 ± 0.13 for the ORD 
PNe, and 3.32 ± 0.09 for the H II regions. The fact that the CRD PN 
with the highest metallicity, NGC 6826, has a value of O/Cl similar 
to the ones found for ORD PNe and H II regions suggests that O 
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Figure 3. Val ues of O/Cl as a function of O/H (left panel) and Cl/H (right panel) for our sample of ORD PNe, CRD PNe, and H 11 regions. The protosolar 
abundances of lLodde^ j2010ll are overplotted with the solar symbol. We use the same symbols as in Figure[T] 


enrichment is no longer important at metallicities close to solar, but 
this should be confirmed by the analysis of other CRD PNe with 
similar metallici ty. Note that N GC 6826 is suspected to have a bi¬ 
nary companion l lMendezll989t) . which might be responsible for its 
carbon enrichment. 

We performed a rough estimation of the amount of oxygen 
produced by low-mass stars and hi gh-mass stars using the yields 
calculated bv IPignatari et a a (l2013h . To do this, we integrated the 
net yields over the Salneter 1 19551) initial mass function in two mass 
ranges: 1.65-5 Mq and 15-60 Mq forZ = 0.02, 1.65-5 Mq and 15- 
25 Mo for Z = 0.01. We obtain that low-mass stars contribute about 
4 per cent and about 10 per cent of the total oxygen produced by a 
single stellar population at Z = 0.02 and 0.01, respectively. These 
percentages are small, but they are likely to increase at lower metal¬ 
licities, and detailed chemical evolution models will be needed to 
find out if low-mass stars are important producers of oxygen in dif¬ 
ferent galaxies. 

Figure 0 offers an explanation for the overabundance 
of oxygen in PNe of the solar neighbourhood found by 
iRodriguez & Delgado-Ingladal ( 1201 ih . This result would be due to 
a combination of two causes: the presence of high metallicity PNe 
in our vicinity for any of the reasons listed above plus oxygen en¬ 
richment in some of the lower metallicity PNe. 


4.5 Nitrogen production in the PNe with CRD and ORD 

Figurel^shows the values of the N/0 abundance ratio as a function 
of the helium abundance. This plot is commonly used to identify 
the descendants of massive progenitors, which are efficient produc¬ 
ers of N and Fie according to the models of AGB evolution. The 
top panel of Figure [5] shows t he results obtained using the ICFs of 
IPelgado-Inglada et al. I ( l2014h . The IF II regions show some disper¬ 
sion in their values of N/O and we noted that it seems to be related 
to their degree of ionization, as measured with 0^^/(0^ + 0^ *^): the 
three H II regions of higher ionization have the largest values of 
N/O. We decided to compare the results with those obtained using 
the classical ICF for nitrogen: N/O = N^/O^. The middle panel of 
Figure 13 shows the abundances implied by this ICF. It can be seen 
that the new N/O ratios of the H II regions have less dispersion and 
a better agreement with the solar value. On the other hand, the val¬ 


ues of N/O in the PNe with higher helium abundances now show a 
tight correlation with the helium abundance. We cannot think of any 
observational effect or bias that would introduce this correlation if 
it were not real, and we conclude that the classical ICF for nitrogen 
seems to be working better, at least for this sample of objects. The 
issue merits more investigation, but since our conclusions do not 
depend on the chosen ICF, we will not discuss this further. 

The bott om panel of Fig ure [3 compares our results with the 
predictions of iKarakal l l201(j|) an mPignatari et alj j2013t) . The se¬ 
quence of triangles with increasing size correspond to models with 
initial metallicities of Z = 0.008 and 0.02 and masses of 1, 1.25, 
1.5, 1.75, 1.9, 2, 2.25, 2.5, 3, 3.5, 4, 4.5, 5, 5.5, 6, an d 6.5 M^, 
(the se quence ends at 6 Mq for Z = 0.008) computed bv iKarakasI 
1 20101). The sequenc e of hexagons correspond to the models of 
IPignatari et alj ( 12013h with initial masses of 1.65, 2, 3, and 5 Mq 
and Z = 0.01 and 0.02. According to the models, most of our 
ORD PNe descend from stars with masses above 4.5 Mq, bom in 
a medium of near-solar composition, whereas the CRD PNe would 
descend from stars with masses below ~4 Mq formed in an envi¬ 
ronment with sub-solar metallicity. Note however, that extra mix¬ 
ing processes, stellar rotation, or the effect of a binary compan¬ 
ion could produce large N/O ratios in PNe with initial masses not 
that high jKarakasiboid : [^sinska et al.ll2013h . In particular, the 
models we are considering imply that all but one of the bulge PNe 
have high-mass progenitors, as discussed by other auth ors (see e.g. 
IWang & Lhil 120071 : iGarcia-Hernandez & Gornvl l2014h. Ho wever, 
this conclusion should be considered with caution. Isuelll ( 12013h 
argues that if the bulge PNe had progenitor masses above 4 Mq, 
we would observe in the bulge bright carbon stars arising from 
stars o f 2-4 Mq, which is not the case. Following the sugges¬ 
tion of iNataf & Gouldl ( 12012h that there is a population of helium- 
enhance d stars in the bulge, similar to those found in some globular 
clusters. iBue il( l2013h constructs a set of AGB models with masses 
between 1.2 and 1.8 Mq, Z ~ 0.02, and pre-enriched in helium, 
finding that they can reproduce the observed abundance ratios. 

Since chlorine is a better proxy for the metallicity than oxy¬ 
gen, we display in Figurej^the results for the nitrogen abundances 
relative to chlorine. We also show t he p redictions from the nu cle- 
osynthesis models of lKarak^ ( l2010l) and IPignatari et al.l ( l2013h . As 
expected, the results for the ORD PNe are almost the same when 
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Figure 5. Values o f N/O as a function of He/H. The protosolai' abun¬ 
dances of lLodde^ ilOlOl) are overplotted with the solar symbol. We use 
the same symbols as in Figure ^ We draw upper limits in the He abun¬ 
dance of four objects where the contribution of neutral helium is signifi- 
ca nt. The top panel shows the ni trogen abundances derived with the ICFs 
of IPelgado-Inglada et al. 1 <20l4 : in the middle panel the nitrogen abun¬ 
dances have been calculated using the classical ICF, N/O = N+/0^. The 
bottom panel s how s the predictions from the nucleosynthesis models by 
iKarakasI ilOlOl) and IPignatari et al] i2013h overplotted for comparison with 
the results of the middle panel. The size of the symbols is larger for higher 
mass progenitors (see the text for more information). 
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Figure 6. Values of N/Cl as a function of He/H. The N abundances have 
been calculate d using the cla ssical ICF, N/O = N^/0^. The protosolar 
abundances of lLodder^ j2010l) are overplotted with the solar symbol. We 
use the same symbols as in Figure]^ We draw upper limits in the He abun¬ 
dance of four objects where the contribution of neutral helium is significant. 


compared to the values in H II regions because their oxygen abun¬ 
dance has not been modified, but the N/Cl values obtained in the 
CRD PNe reveal that most of them are also significantly enriched 
in nitrogen. 


4.6 The C/O abundance ratio 

Figure |7] shows the C/O values as a function of Fle/Fl, both de¬ 
rived from RLs. For three of the ORD we do not have an esti¬ 
mate of the C/O abundance ratio and they do not appear in this 
figure. In principle, RLs and CELs m ight be tracing different nebu¬ 
lar phases (see e.g. lLiu et al.ll20045) . but the valu es of C/O derived 
using RLs and CELs are similar in most cases jLiu et al.ll200^: 
IPelgado-Inglada & Rodrfguezll2014l) . 

We can see that the value derived for C/O does not provide 
a direct indication of the type of dust grains present in a PN. This 
can be attributed to the uncertainties involved in its determination, 
since the C/O value is based on weak RLs and requires an ICF. In 
fact, the dispersion of our Monte Carlo error simulations shows that 
H II regions and ORD PNe have C/O values compatible with the 
expectation that C/O < 1 in these objects and the values of CRD 
PNe are also compatible with C/O > 1, as expected from their dust 
features. The average value of C/O for CRD PNe is also higher 
than the one shown by the other objects, which cluster around the 
solar value. Another effect that could affect the observed value of 
C/O is depletion of carbon or oxygen into dust grains. For these 
reasons, the dust features observed in the infrared spectra are likely 
to provide more reliable information on whether the value of C/O 
is below or above 1 in each object. 

If we now examine the values of He/H (including those listed 
in Table [3 for the three ORD PNe that have no estimate of C/O), 
we can see that excluding the low-metallicity object from the halo, 
the ORD PNe have helium abundances which are systematically 
higher than those found for CRD PNe and H II regions. This is 
consistent with our conclusions above that they had massive (or 
pre-enriched) stellar progenitors, since AGB models predict that 
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Figure 8. Values of C/0 as a function of O/H derived from RL s for the 
studie d PNe and the H II regions. The protosolar abundances of lLodder^ 
<2O10l) are oveiplotted with the solar symbol. We use the same symbols as 
in Figure[2 Note that for three of the ORD PNe we could not estimate the 
value of C/O and these objects are not plotted in this figure. 


helium enrichment takes place for the stars with the highest masses. 
This is illustrated in the right panel of the figure, where we plot the 
same models we used for Figure]^ 

In Figure[8l we also show the derived values of C/O but now as 
a function of the oxygen abundances implied by RLs (listed in Ta- 
bleH). Since the correct values for the nebular oxygen abundances 
are likely to be intermediate between the ones implied by CELs 
and RLs (or close to either one), the O/H values derived from RLs 
complete our picture of the distribution in metallicity of the sam¬ 
ple objects and their comparison with the Sun. The O/H values for 
H II regions and ORD PNe should be corrected upward by 0.10- 
0.15 dex to take into account the presence of oxygen in silicates 
and o xides within dust grains ( IWhittelil20ld : |Peimbert & Peimberli 
l2010li . With this correction, the oxygen abundance in the H II re¬ 
gions would still be below solar for the results based on CELs, but 
would agree both with the solar abundance and with the abundances 
found for nearby B stars jNieva & Przvbillj|2012h for the results 
based on RLs. On the other hand, the oxygen abundances derived 
with CELs imply that the ORD PNe of the solar neighbourhood 
have a chemical composition which is close to solar, whereas RLs 
imply that these objects are significantly more metal-rich than the 
Sun. Both RLs and CELs imply that several nearby PNe are more 
metal-rich than the H II regions, a result confirmed by the derived 
values of Cl/H. 

Eigure [8] a lso shows the pred i ctions of the AGB models of 
iKarakaj l l20ldl and IPignatari et all ( 1201 3h . In both sets of models 
the results for the lower mass objects indicate the initia l value s of 
O/H. As we commented above, the models of lKarak^ | |201(]|) for 
the most massive stars with Z = 0 .008 predict significan t oxygen 
destruction, whereas the models of IPignatari et alj l l2013ll with ini¬ 
tial masses of 2-3 Mo are enriched both in carbon and oxygen. As 
discuss ed above, our results fo r the CRD PNe agree with the predic¬ 
tions of IPignatari et alj ( l2013h . at least for Z = 0.01. An extended 
analysis of other CRD PNe might provide more information. 


5 CONCLUSIONS 

We have studied the chemical composition of 20 Galactic PNe with 
high quality infrared and optical spectra. The sample includes one 
halo PN, eight bulge PNe, and eleven disc PNe. We have also se¬ 
lected from the literature seven H II regions from the solar neigh¬ 
bourhood with the best available optical spectra.Using the infrared 
dust features we have classified the PNe in two groups: PNe with 
ORD and PNe with CRD. The first ones show amorphous and/or 
crystalline silicates in their spectra whereas the second ones show 
the broad features at 11 and/or 30 jiva associated with SiC and MgS. 

We have computed the physical conditions and the chemi¬ 
cal abundances following the same procedure in all the objects. 
The abundances derived here are based on the best av ailable ICEs, 
mainly those derived bv IPelgado-Inglada et al. I ( l2014h . and the un¬ 
certainties in the element abundances include the uncertainties in¬ 
troduced by the physical conditions, the ionic abundances, and the 
ICEs. 

We find a tight correlation between the abundances of Cl and 
Ar, in agreement with a lockstep evolution of both elements. The 
metallicity in our sample, traced by Cl/H, covers one order of 
magnitude, and there are significant differences in the Cl/H ratio 
among the three groups of objects studied here. In general, the PNe 
with CRD have low metallicity, the H II regions have intermedi¬ 
ate metallicities, and the PNe with ORD cover the whole range of 
metallicity (with some of them showing the highest abundances, 
even higher than the solar ones). This result indicates that the pro¬ 
genitors of the PNe with CRD were formed in a sub-solar metallic¬ 
ity medium (with the exception of NGC 6826, with a value of Cl/H 
similar to the H II regions) whereas most of the PNe with ORD 
were formed in a near-solar metallicity medium (with the excep¬ 
tion of the halo PN DdDm 1). 

The neon abundances at a given metallicity show a relatively 
large dispersion. The dispersion could be due to uncertainties in the 
ICE but we cannot rule out that some of the PNe are neon enriched. 

The oxygen abundances do not always trace the metallicities 
of the objects. The 0/Cl values in the Sun, the PNe with ORD, 
and the H II regions are similar, but most of the PNe with CRD 
show 0/Cl values that are higher by a factor of two. The exception 
is NGC 6826, the CRD PN with higher metallicity in our sample. 
As we mentioned above, the progenitors of most of the PNe with 
ORD were formed in an environment with near-solar metallicity, 
where the nucleosynthesis process that can change the initial oxy¬ 
gen abundance are less efficient. This is in agreement with the sim¬ 
ilar values of O/Cl found in these PNe and in the H II regions. The 
low Cl/H abundances and the high 0/Cl values found in all but one 
of the PNe with CRD provide evidence of oxygen production in the 
progenitor stars of PNe with CRD. Non-stand ard nucleosynthesis 
models, such as those bv IPignatari et al.l(l2013h . that include an ex¬ 
tra mixing mechanism, predict a significant production of oxygen 
via the third dredge-up in stars with masses around ~ 2-3 Mq, ex¬ 
plaining our results. The relevance for galactic chemical evolution 
of oxygen production in low mass stars deserves further study using 
chemical evolution models. 

The comparison of the He/H, C/O, and N/0 abundance ra¬ 
tios derived for the PNe with the predictions of nucleosynthesis 
models, suggests that the PNe with CRD descend from stars with 
masses in the range 1.5-3 Mo whereas the PNe with ORD and the 
highest abundances of nitrogen and helium would arise from stars 
with higher initial masses. The theoretical models indicate that they 
descend from stars with masses > 4.5 Mo, but this limit would be 
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Figure 7. Values of C/0 as a function of He/H for the studied PNe and the H II regions. The protosolar abundances of lLodde^ jlOldi are overplotted with the 
solar symbol. We use the same symbols as in Figure[T] Note that for three of the ORD PNe we could not estimate the value of C/0 and these objects are not 


plotted in this figure. 


smaller with models that include stellar rotation, extra mixing pro¬ 
cesses or pre-enrichment. 

Our results indicate that Ar and Cl are the best metallicity in¬ 
dicators for PNe. At any metallicity their abundances reflect the 
composition of the interstellar medium where the PN progenitors 
were formed. Moreover, the ICFs adopted here for Ar and Cl can 
be used when only CP^ and Ar^^ lines are observed, making easier 
a homogeneous calculation of PN metallicities. High resolution and 
deep spectra are needed, especially of extragalactic PNe, to delve 
in this issue and calculate, for example, metallicity gradients using 
Cl or Ar. Besides, the obtention of high quality optical spectra in 
which we can measure Cl and Ar lines will provide more reliable 
estimates of the electron densities. The observed extragalactic PNe 
are often the brightest ones, with high electron densities, where the 
commonly used density diagnostic ratio based on [S II] lines is not 
sensitive. 


ACKNOWLEDGEMENTS 

We are grateful to the anonymous referee for valuable comments 
and suggestions. The authors thank L. Carigi, L. Hernandez- 
Martmez, and M. Pignatari for fruitful discussions and comments. 
G.D.-I. gratefully acknowledges a DGAPA postdoctoral grant 
from the Universidad Nacional Autonoma de Mexico (UNAM). 
C.M. and G.S. acknowledge support from the following Mexi¬ 
can projects: CB-2010/153985, PAPIIT-IN105511, and PAPIIT- 
INl 12911. M.P. is grateful for the finantial support provided by 
CONACyT grant 129753. M.R. and G.D.-I. acknowledge support 
from Mexican CONACyT grant CB-2009-01/131610. This work 
has made use of NASA’s Astrophysics Data System, and the SIM- 
BAD database operated at CDS, Strasbourg, France. 

REFERENCES 

Allende Prieto C., Lambert D. L., Asplund M., 2002, ApJ, 573, 
L137 


Asplund M., Grevesse N., Sauval A. J., Scott P, 2009, ARA&A, 
47, 481 

Barker T., Cudworth K. M., 1984, ApJ, 278, 610 
Bernard-Salas J., Peelers E., Sloan G. C., Gutenkunst S., Matsuura 
M., Tielens A. G. G. M., Zijlstra A. A., Houck J. R., 2009, ApJ, 
699, 1541 

Buell J. F., 2013, MNRAS, 428, 2577 

Ciardullo R., Bond H. E., Sipior M. S., Fullton L. K., Zhang C.-Y, 
Schaefer K. G., 1999, AJ, 118, 488 
Clegg R. E. S., Peimbert M., Torres-Peimbert S., 1987, MNRAS, 
224, 761 

Davey A. R., Storey P. J., Kisielius R., 2000, A&AS, 142, 85 
Delgado-Inglada G., Morisset C., Stasihska G., 2014, MNRAS, 
440, 536 

Delgado-Inglada G., Rodriguez M., 2014, ApJ, 784, 173 
Depew K., Parker Q. A., Miszalski B., De Marco O., Erew D. J., 
Acker A., Kovacevic A. V., Sharp R. G., 2011, MNRAS, 414, 
2812 

Esteban C., Garcla-Rojas J., Peimbert M., Peimbert A., Ruiz 
M. T., Rodriguez M., Carigi L., 2005, ApJ, 618, L95 
Esteban C., Peimbert M., Garcla-Rojas J., Ruiz M. T., Peimbert 
A., Rodriguez M., 2004, MNRAS, 355, 229 
Galavis M. E., Mendoza C., Zeippen C. J., 1997, A&AS, 123, 159 
Garcia-Hernandez D. A., Gorny S. K., 2014, A&A, 567, A12 
Garcla-Rojas J., Esteban C., Peimbert A., Rodriguez M., Peimbert 
M., Ruiz M. T, 2007, Rev. Mex. Astron. Astrohs., 43, 3 
Garcla-Rojas J., Esteban C., Peimbert M., Rodriguez M., Ruiz 
M. T., Peimbert A., 2004, ApJS, 153, 501 
Garcla-Rojas J., Esteban C., Peimbert M., Costado M. T., Ro¬ 
driguez M., Peimbert A., Ruiz M. T., 2006, MNRAS, 368, 253 
Garcla-Rojas J., PenaM., Morisset C., Delgado-Inglada G., Mesa- 
Delgado A., Ruiz M. T., 2013, A&A, 558, A122 
Gorny S. K., Chiappini C., Stasihska G., Cuisinier F, 2009, A&A, 
500, 1089 

Gorny S. K., Stasihska G., Escudero A. V., Costa R. D. D., 2004, 
A&A, 427, 231 

Henry R. B. C., Kwitter K. B., Balick B., 2004, AJ, 127, 2284 
Henry R. B. C., Kwitter K. B., Dufour R. J., Skinner J. N., 2008, 


















14 Delgado-Inglada et al. 


ApJ, 680, 1162 

Herwig R, Bloecker T., Schoenberner D., El Eid, M., 1997, A&A, 
324, L81 

Izotov Y. 1., Stasinska G., Meynet G., Guseva N. G., Thuan T. X., 
2006, A&A, 448, 955 
Jenkins E. B., 2009, ApJ, 700, 1299 
Karakas A. 1., 2010, MNRAS, 403, 1413 
Karakas A. 1., Lattanzio J. C., 2003, RASA, 20, 393 
Karakas A. 1., Lattanzio J. C., 2014, preprint (arXiv: 1405.0062) 
Kingsburgh R. L., Barlow M. J., 1994, MNRAS, 271, 257 
Kisielius R., Storey R J., Ferland G. J., Keenan F. R, 2009, MN¬ 
RAS, 397, 903 

Krueger T. K. Czyzak S. J., 1970, Rroc. R. Soc. Lond. A, 318, 531 
Leisy R, Dennefeld M., 2006, A&A, 456, 451 
Liu Y, Liu X.-W., Barlow M. J., Luo S.-G., 2004b, MNRAS, 353, 
1251 

Liu Y, Liu X.-W., Luo S.-G., Barlow M. J., 2004a, MNRAS, 353, 
1231 

Liu X.-W., Luo S.-G., Barlow M. J., Danziger 1. J., Storey R. J., 
2001, MNRAS, 327, 141 

Liu X.-W., Storey R. J., Barlow M. J., Danziger 1. J., Coben M., 
Bryce M., 2000, MNRAS, 312, 585 
Lodders K., 2010, in Rrinciples and Rerspectives in Cosmocbem- 
istry, eds. A. Goswami & B. E. Reddy, Astropbys. Space Sci. 
Rroc. (Springer-Verlag), 379 

Luridiana V., Morisset C., Sbaw R. A., 2014, preprint 
I arXiv: 1410.66621 
Marigo R, 2001, A&A, 370, 194 

McLaughlin B. M., Bell K. L., 2000, Journal of Rbysics B Atomic 
Molecular Rbysics, 33, 597 
Mendez R. H., 1989, lAUS, 131, 261 
Mendoza C., Zeippen C. J., 1982, MNRAS, 198, 127 
Milingo J. B., Kwitter K. B., Henry R. B. C., Souza S. R, 2010, 
ApJ, 711,619 

Munoz Burgos, J. M., Locb, S. D., Ballance, C. R, Boivin, R. F, 
2009, A&A, 500, 1253 
Nataf D. M., Gould A. R, 2012, ApJ, 751, 39 
Nieva M.-F, Rrzybilla N., 2012, A&A, 539, 143 
Reimbert A., Reimbert M., 2010, ApJ, 724, 791 
Reimbert A., Reimbert M., Ruiz M. T., 2005, ApJ, 634, 1056 
Requignot D., Walsb J. R., Zijlstra A. A., Dudziak G., 2000, A&A, 
361, 1 

Rignatari M., et al., 2013, preprint (arXiv: 1307:6961) 
Rodobedova L.L, Kelleber D.E., Wiese W.L., 2009, JRCRD, 38, 
171 

Rorter R. L., Ferland G. J., Storey R. J., Detiscb M. J., 2012, MN¬ 
RAS, 425, L28 

Rorter R. L., Ferland G. J., Storey R. J., Detiscb M. J., 2013, MN¬ 
RAS, 433, L89 

Quireza C., Rocba-Rinto H. J., Maciel W. J., 2007, A&A, 475,217 
Ramsbottom C. A., Bell K. L., Keenan F. R, 1997, MNRAS, 284, 
754 

Rodriguez M., Delgado-Inglada G., 2011, ApJ, 733, 50 
Rubin, R. H. 1986, ApJ, 309, 334 
SalpeterE. E., 1955, ApJ, 121, 161 

Sbarpee B., Williams R., Baldwin J. A., van Hoof R. A. M., 2003, 
ApJS, 149, 157 

Simon-Dlaz S., Stasinska G., 2008, MNRAS, 389, 1009 
Stasinska G., Rena M., Bresolin F, Tsamis Y. G., 2013, A&A, 
552, 12 

Stasinska G., Tylenda R., Acker A., Stenbolm B., 1991, A&A, 
247, 173 


Stasinska G., et al., 2010, A&A, 511, AA44 
Storey R. J. 1994, A&A, 282, 999 
Storey R. J., Hummer D. G., 1995, MNRAS, 272, 41 
Storey R. J., Zeippen C. J., 2000, MNRAS, 312, 813 
Storey R. J., Socbi T., Badnell N. R., 2014, MNRAS, 441, 3028 
Tayal S. S., 2011, ApJS, 195, 12 
Tayal S.S., Gupta G.R., 1999 ApJ 526, 544 
Tayal S. S., Zatsarinny O., 2010, ApJS, 188, 32 
Tylenda R., Acker A., Stenbolm B., 1993, A&AS, 102, 595 
Tsamis Y. G., Barlow M. J., Liu X.-W., Danziger 1. J., Storey R. J., 
2003, MNRAS, 345, 186 
Wang W., Liu X.-W., 2007, MNRAS, 381, 669 
Waters L. B. F. M., et al., 1998, A&A, 331, L61 
Wesson R., Liu X.-W., Barlow M. J., 2005, MNRAS, 362, 424 
Wbittet D. C. B., 2003, Dust in Galactic Environment. loR Rub- 
lisbing, Bristol 

Wbittet D. C. B., 2010, ApJ, 710, 1009 
Zeippen C. J., 1982, MNRAS, 198, 111 
Zbang C. Y, 1995, ApJS, 98, 659 


